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Abstract 

The two-dimensional unsteady coupled Burgers' equations with moderate to severe 
gradients, are solved numerically using higher-order accurate finite difference schemes; 
namely the fourth-order accurate compact ADI scheme, and the fourth-order accurate 
Du Fort Frankel scheme. The question of numerical stability and convergence are 
presented. Comparisons are made between the present schemes in terms of accuracy 
and computational efficiency for solving problems with severe internal and boundary 
gradients. The present study shows that the fourth-order compact ADI scheme is 
stable and efficient. 

1 Introduction 

The basic flow equations describing unsteady transport problems form a parabolic hyper- 
bolic system of partial differential equations. The interaction of the nonlinear convective 
terms and the dissipative viscous (or dispersion) terms in these equations can result in 
relatively severe gradients in the solution. Also, the accuracy of the numerical solution and 
the computational efficiency are highly dependent on the numerical methods used to solve 
this kind of partial differential equations. Standard three-point finite difference methods 
of approximating spatial derivatives may work well for smooth solutions, but they fail 
when severe gradients or discontinuities are present, which are common in the shock wave 
problems [1, 2, 3]. Lower-order accurate finite difference methods, such as upwinding- 
type finite differences, can be a remedy for the numerical oscillations and dispersions. 
However, they have a large amount of "numerical viscosity" that smooths the solution 
in much the same way that physical viscosity would, but to an extent that is unrealistic 
by several orders of magnitude, see Sharif and Busnaina [4]. Standard four-point finite 
difference methods, such as Leonard's method [5], are good in their higher-order accura- 
cies and in reducing numerical smearing effects. But, they are plagued by their generation 
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of spurious oscillations or overshoots in the neighborhood of discontinuities and lack ac- 
curacy, as reported by Liu et al. [6], and the present author [7]. TVD finite difference 
schemes [8, 9], guarantee oscillation- free solutions, but they are limited to second-order 
accuracy. Third-order accurate TVD schemes are reported by Gupta et al. [10], and by 
Liu et al. [6]. However, in their studies, nonlinearities were either not present or played 
a minor role. Therefore, the perfect numerical methods should possess both higher-order 
accuracy and sharp resolution of discontinuities without excessive smearing. Moreover, 
higher-order accurate numerical methods are attractive for problems with long compu- 
tational time or with required higher accuracy solutions. As mentioned by Orszag [11], 
that at the cost of slight additional computational complexity, the fourth-order schemes 
achieve results in the 5% accuracy range with approximately half the spatial resolution 
in each space direction compared with the second-order schemes (i.e. a factor 8 fewer 
grid points in three dimensions). But, the objection to the standard higher-order schemes 
comes from the additional nodes necessary to achieve the higher-order accuracy. This 
precludes the use of implicit methods since the obtained matrix is not of tridiagonal form, 
and it is necessary to use fictitious nodes for the boundary conditions. Also, they do not 
allow easily for non-uniform grids, unless at the expense of the order of accuracy. On the 
other hand, the compact schemes that treat the function and its necessary derivatives as 
unknowns at the grid nodes, like the Pade scheme [40], are fourth-order accurate, and 
compact in sence that they reduce to tridiagonal form. The compact schemes generally 
consist of finite difference schemes which involve two or three grid points. The three-point 
schemes fall into two classes. The first class consists of methods which are fourth-order 
accurate for uniform grids, such as Kreiss scheme [11, 12], the Mehrstellen method [13], the 
operator compact implicit scheme [14-16], and the Hermitian finite difference method of 
Peters [17]. The second class consists of methods that allow variable grids such as the cubic 
spline methods of Rubin and Graves [18-21], and the Hermitian finite difference method 
of Adam [22, 23]. For incompressible viscous flow problems, there has been several work 
on the construction of compact schemes for the incompressible Navier-Stokes equations, 
see for example, [24-33]. The most noted ones include the work of Gupta [24, 25]. He 
introduced a compact fourth-order finite difference scheme with three nodal points for the 
convection diffusion equations. His scheme does not seem to suffer excessively from spu- 
rious oscillatory behavior or numerical viscosity, and it converges with standard methods 
such as Gauss Seidel or SOR regardless of the diffusion [34, 35]. However, Yavneh [36], 
in his analysis of Gupta's scheme for the 2-D convection-diffusion equation with Dirichlet 
boundary conditions, has concluded that Gupta's compact scheme does not suffer from a 
cross-stream artificial viscosity, but it does include a streamwise artificial viscosity that is 
inversely proportional to the natural viscosity. For unsteady viscous flow cases, Weinan 
and Liu [37] introduced an modified version of Gupta's compact scheme based on the 
vorticity-stream function formulation. The requirement of compactness of this scheme is 
slighly relaxed (for the convection terms) so that the resulting scheme is simple, easy to 
implement, and with small phase error. They were successful to compute the driven cavity 
problem at high Reynold numbers upto 10 6 on 1024 2 grid, and they have concluded that 
the fourth-order compact methods are comparable in accuracy with the spectral methods 
for most problems of practical interests. For turbulent fluid flow, where there is a range of 
space and time scales, Lele [38] had introduced a series of higher-order compact schemes 
that are generalization of the Pade scheme with three nodal points, and with an improved 
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representation of the shorter length scales. He applied them to the evolution of supersonic 
shear layers. 

The disadvantage of the above higher-order compact schemes involving three nodal 
points is that the boundary conditions are no longer sufficient and they do not allow 
easily for non-uniform grids, unless at the expense of the order of accuracy. Another 
disadvantage of some compact schemes is the complexity of the resulting nonlinear finite 
difference equations and the associated difficulty in solving them efficiently. On the other 
hand, the compact scheme with two nodal points, like second-diagonal Pade scheme, is 
fourth-order accurate even for non-uniform spatial grids, and no fictitious points neither 
extra formula are needed for Dirichlet boundary conditions, as discussed by White [39] 
and Keller [40]. Also, Liniger and Willoughby [41] studied the numerical solutions of 
stiff systems of ordinary differential equations, that are encountered in many areas of 
applied mathematics, using compact two-point implicit methods. They introduced three 
main compact schemes with different order of accuracy, and with some very favorable 
properties. In particular, their schemes have A-stability in the sense of Dahlquist [55] 
and they account for the exponential character of the rapidly decaying solutions directly, 
which are referred as exponential fitting methods. 

Inspite of many articles have appeared in the literature concerning the applications 
of the higher-order accurate schemes including the compact schemes to fluid dynamics 
problems, there is no much works done in the area of application of two-point compact 
schemes [39-44], like the fourth-order accurate second-diagonal Pade approximation, to 
multi-dimensional cases. This is the main objective of present study, where we study 
the feasibility of extending the two-point compact scheme to solve the unsteady two- 
dimensional coupled Burgers' equations. They take the following form: 

du du du d 2 u d 2 u _ 



dt 




(1.1) 



dv 
dt 




(1.2) 



Xq < X < XjV, 



yo < y < vm, 



t > 



with the initial conditions 



u(x,y,0) = u (x,y), 



v(x,y,0) = v (x,y) 



(1.3) 



and the Dirichlet boundary conditions 



u{x ,y,t) = ui(y), 
u{x,y ,t) = u 3 {x), 



u{xN,y,t) = u 2 {y), 
u(x,y M ,t) = u 4 (x), 



(1.4) 



v{x ,y,t) = ui(y), 
v(x,y ,t) = v 3 (x), 



v(x N ,y,t) = v 2 (y), 
v(x,y M ,t) = vt(x), 



(1.5) 



where v is equal to 1/Re > 0, and Re is the Reynolds number. Coupled Burgers' equations 
are used to model many practical transport problems, such as vorticity transport, hydrody- 
namic turbulence, shock wave theory, wave processes in thermoelastic medium, transport 
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and dispersion of pollutants in rivers, and sediment transport, see references [45-47]. The 
coupled Burgers' equations are an appropiate form of the Navier-Stokes equations. They 
have the same convective and diffusion form as the incompressible Navier-Stokes equa- 
tions. Fletcher [48] made a comparison of finite element and finite difference methods 
with different orders of accuracy for solving the two-dimensional Burgers' coupled equa- 
tions, and he concluded that the five-point finite difference scheme is the most efficient 
scheme. Boonkkamp and Verwer [49] have used the extrapolated Odd-Even Hopscotch 
scheme for solving the inhomogeneous two-dimensional coupled Burgers' equations, but 
their solution exhibited wiggles. Arminjon and Beauchamp [50] have concluded that the 
finite element method is efficient compared to the other methods, namely the method of 
lines and Runge-Kutta-type method in solving the above Burgers' equations. Jain and 
Raja [51] have used splitting-up technique to reduce the problem to a sequence of tridiag- 
onal systems. Jain and Lohar [52] used spline locally one dimensional (SOLD) algorithm 
for solving coupled nonlinear parabolic equations. EL-Zoheiry and EL-Naggar [53] used 
the spline alternating direction implicit (SADI) method for solving the two-dimensional 
Burgers' equations. EL-Naggar [54] has presented a mixed implicit-explicit two levels algo- 
rithm that is based on SADI method for solving the coupled Burgers' equations. However, 
his results exhibited oscillations. 

In the present study, higher-order accurate two-point compact alternating direction 
implicit algorithm (CADI) is introduced to solve the two-dimensional unsteady coupled 
Burgers' equations, for problems with moderate to severe internal and boundary gradients. 
The algorithm has the following features: 

(1) it results in finite difference equations that involve only two-nodal points and there- 
fore is formally fourth-order accurate on all grid points, even for non- uniform grids; 

(2) it has ^4-stability in the sense of Dahlquist, and accounts for the exponential character 
of rapidly varying solutions [41]; 

(3) it utilizes Newton's method for linearization with a quardratic convergence; 

(4) it requires only the given Dirichlet boundary conditions; 

(5) the algorithm is simple and easy to implement. 

In short, the present method is the natural extension of A-stable fourth-order accurate 
second-diagonal Pade approximation to solve multi-dimensional flow problems with mod- 
erate to severe gradients. Comparison of the present scheme with the fourth-order Du 
Fort Frankel scheme is made in terms of accuracy and computational efficiency, which 
show that the fourth-order compact ADI scheme is stable and efficient. 

2 The numerical schemes 

In this section, the present numerical schemes, namely the fourth-order accurate two-point 
compact ADI scheme and the fourth order Du Fort Frankel scheme, are derived for the 
two-dimensional unsteady coupled Burgers' equations (1). 



On the Fourth-Order Accurate Compact ADI Scheme 



17 



2.1 Fourth-order accurate Du Fort Frankel scheme 

Let the interval [xq, xn] be discretized into N grid steps of size Ax, where Ax = (xj — Xj_i), 
i is an index of any grid-point in x direction. Similarly, the interval [yo,UM] is discretized 
into M grid steps of size Ay, where Ay = (yj — yj-i), j is an index of any grid point in 
y-direction, and n is an index for the temporal grid point. The explicit form of the Du 
Fort Frankel scheme for the two-dimensional coupled Burgers' equations (1), using Kreiss 
fourth-order accurate approximations [11] for the spatial derivatives, takes the following 
form: 



n+l 



n-l 



2At 



-u%D x 



1 Ax V 
!--6-^ 



Ay% 2 



Ar 2 
12 x 



(2.1) 



where 



ip = [u, v] 



(2.2) 



(2.3) 



(2.4) 



Define c™ 



u"^ c™ 
« Ax' y 



(L 



v— — ^, and cL 
Ax 2 y 



At 
Ay 2 



At 

v%a— — , to be the local Courant numbers in x and y directions, 
%3 Ay 



The above equation represents the fourth-order accurate 



leap-frog scheme for equations (1), and in order to obtain the final form of the fourth- 
order accurate explicit Du Fort Frankel scheme for the 2-D unsteady coupled Burgers' 
equations (1), the center node value (ipij) in the diffusion terms in equation (2.1) are 
replaced by their average at time-levels (n — 1) and (n + l), giving: 



(2.5) 



where 



A = (1 -2.5d x -2.5d y )/Q, 
C = (-8c™ + 16d x )/6Q, 
E = -(<% + d x )/6Q, 
G = (-8c™ + Wd y )/6Q, 
L = -(c^ + d y )/6Q, 



B = (c x - d x )/6Q, 
D = (8c™ + 16d x )/6Q, 
F = (c™ - dy)/6Q, 
H = (8c™ + 16d y )/6Q, 
Q = (l + 2.5d x + 2.5dy). 



(2.6) 
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2.2 Fourth-order accurate compact ADI scheme 



Liniger et al. [41] have introduced the following linear one step formulas for <p(x) containing 
real free parameters (a Sz b): 

Ax Ax 2 
[<t>i+i -4>i\-~2 K 1 + + (!-«)&*] + i( b + a ) ( l ) xx i+1 ~ (b- a)(f) XXi ] = e T ,(2.7) 



Ax 3 f 1 



J [2 f - 2(1 - a)i + (b- a)] ^ 0(x + CAx)de 



(2.8) 



For the case of a > and 6 = 1/3, the resulting formula has a third-order accuracy. 
Moreover, for the case of a = and 6 = 1/3, the formula has a fourth-order accuracy, 
which is known as the two-point second-diagonal Pade approximation: 

At Ax 2 
b i+1 - <t>i\ - — [4>x l+1 + <t> Xi \ + -T7r[®xxi+i ~ ^xxi] = er, (2.9) 



Ax ; 
~24 



\ ~ —[4>x t+1 + <f> Xi ] + -J^[$xx i+1 - 
rl Q5 

J e 2 «-l) 2 ^5^ + eAxK, 



(2.10) 



where <j) x , 4> xx are the first and the second derivatives of the function 4>(x). Using the 
above scheme and an ADI-type time marching procedure for the temporal derivative, 
the compact alternating direction implicit algorithm (CADI) for the coupled Burgers' 
equations (1), are obtained by first rewriting these equations as follows: 



u 




[yu x — 0.5n 2 )a; 


V 


t 


{yv x - uv) x + vu x 



+ 



{vUy — UV)y + UVy 

(uv y — 0.5w 2 ) y 



(2.11) 



then the ADI - type time marching procedure requires, in one full time step, the solution 
of: 

x- sweep 

n+0.5 

+ 



v t 



n+0.5 



[yu x — 0.5m ) x 
(vv x - uv) x + vu x 



9i 

92 



n+0.5 



F 
G 



n+0.5 



y-sweep 



u„ 



J y J 



n+l 



n+1 



'\VUy — UV)y + UVy 

(yvy — 0.5f 2 ) y 



n+l 



+ 



fl 
h 



n+0.5 



H 
T 



n n+l 



(2.12a) 
(2.12b) 

(2.13a) 
(2.13b) 



where, \fi, , [di,92] T are the components of the first and the second matrices in the 
right hand side of equation (2.11). The solution procedure consists of solving, first, equa- 
tions (2.12) in the solution vector [U, V, F, G] at time level n+0.5, (x-sweep), then solving 
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equations (2.13) in the solution vector [U, V, H, T] T at time level n + 1, (y-sweep). Noting 
that a = 1 /At, and in order to apply the compact scheme to the solution in the x-sweep, 
a vector Q and its derivatives with respect to x, for Burgers equations (2.12), have been 
defined as follows: 



/«i+0.5 



vu x — 0.5ir 
vv x — uv 
vu 
uv 



n+0.5 



uF - 0.5C/ 2 
vG-UV 
vU 
vV 



n+0.5 



(2.14) 



n+0.5 



Q n+ - 





n+0.5 


~9i ' 
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' aU 


n+0.5 


gi + aU 
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v t - vu x 




92 




aV -VF 




g 2 + aV 




vF 









vF 









uG 


ij 


_ 


ij 


vG 


ij 





ij 



, (2.15) 



Q 



n+0.5 
xxij 



F t 

G t - vF x - Fv x 
u t + uu x 
v t + uv x 



-i n+0.5 




n 




~9l x ' 






92 x 






91 




J ij 


.92 _ 


ij 



aF 


n+0.5 


~aF + g lx ' 


n 


aG - VF X - FG 




aG + g 2x 




aU + UF 




aU + gi 




aV + UG 


ij 


_ aV + g 2 


ij 



(2.16) 



where gi, g 2 , gi x , g 2x and F x in equations (2.15)-(2.16) are approximated by fourth-order 
accurate finite differences. Having substituted the vector Q and its derivatives into the 
above two-point second-diagonal Pade approximation, equation (2.9), by replacing fa by 
the vector Q, we have four non-linear coupled finite difference equations in the solution 
vector [U, V, F, G] T . Newton's method is used to linearize the equations, and the numer- 
ical solution is obtained by iteration. The resulting linearized equations form a block 
tridiagonal matrix system of order N, as in the following form: 



OiSi-i + bi5i + CiS i+ i = ?i 



1,2, 



(2.17) 



where a^, bi, and q are block matrices of order four, 5 = [8U, SV, 8F, 5G] T is the change 
in the solution vector, and f is the right handside vector, each of order four. At each 
iteration, the LU-factorization algorithm is used to obtain the solution of the system (2.17). 
Similarly, the solution procedure of the Burgers' equation (2.11) in the y-sweep, using 
equations (2.13). 



2.3 Numerical stability limits 

The implicit formulation of the two-point compact scheme to the Burgers' equations is 
always unconditionally stable. In this case, the accuracy of the numerical solution depends 
on the size of the discretizations, and higher accuracy can be obtained by finer discretiza- 
tion. Moreover, the present higher-order scheme allows us to use large discretization in 
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comparison with the second-order schemes. Also, it is well known that, for the convection 
diffusion equation, the leap-frog scheme is unconditionally unstable, while the Du Fort 
Frankel scheme has a stability limit (c < 1 & d > 0) [56]. Therefore, it is necessary 
to use Von Neumann stability analysis to define the stability limit of the fourth-order 
accurate Du Fort Frankel scheme to Burgers' equations. In the present study, it will be 
sufficient to examine the stability of the fourth-order Du Fort Frankel scheme for one 
of the above coupled Burgers' equations (2.5), say for u-component. Let the numerical 
solution U(x, y, t) be represented by a finite Fourier series, and for linear stability, we can 
examine the behaviour of a single term of the series, as follows: 



u(iAx,jAy,nAt) = G(nAt)e^ kxiAx+kyjAy ^ , 



(2.18) 



where G(t) is the amplitude function at time-level n of this term whose wave numbers in 
the x and y directions are k x and k y , and / = >/— T. Defining the x and y phase angles as 
&x = k x Ax and 9 y = k y Ay, then, equation (2.18) becomes: 



Qn e I[ie x +jO y ] 



(2.19) 



Substituting (2.19) into the first equation of (2.5), we obtain a quadratic equation for the 
amplification factor (, its solution is: 



c 



Qn+1 

G n 



A ± VA 2 + 4A 



(2.20) 



and 



A = — | [(16 cos 6 X — cos 29 x )d x + (16 cos y — cos 20 y )d y \ 
+/[(-8 sin 9 X + sin 20 x )c x + (-8 sin 9 y + sin 26 y )c y ] } , 



(2.21) 



where A and Q are defined by equations (2.6). For the special case of d x = d y = d and 
c x = c y = c, the modulus of the amplification factor £, defined by the following equation: 



X(c,d,e x ,9 y ) = max 



A + VA 2 + AA 



A - \/A 2 + AA 



(2.22) 



has been computed for different values of |c| and d and plotted, as shown in Fig. 1. This 
shows that the fourth-order accurate Du Fort Frankel scheme is unstable for the range 
(0.35 < |c| < 1.0). E.g. x(l,0.5,7r/2,7r/2) = 1.77, x (l, 0.5, vr, vr) = 1.29, x (0.5, 0.5, tt/2, 
7r/2) = 1.14, and x(0-5, 0.5, it, it) = 1.29. For small values of |c| and d (|c| < 0.35), the 
instability only occurs for phase angles 6 X and 8 y close to it. Moreover, for smaller values 
of d (d < 0.1), the scheme has a neutral stability (\ = 1.0). E.g. x(0.25, 0.5, n/2, tt/2) = 
0.87, x(0.25,0.5,7r,7r) = 1.29, x(0-25, 0.01, vr, vr) = 1.0. Concerning the consistency of 
the present schemes, both of the schemes are consistent with the original differential 
equation (1). The finite difference equation using the compact ADI scheme is consistent 
in sense that the local truncation error, er = O [Ax 5 , At Ax 2 , At Ax] tends to zero as 
At and Ax tend to zero. For Du Fort Frankel scheme equation (2.5) whose truncation 
error, ex = 0[At 2 , (At/ Ax) 2 , Ax 4 ], the consistency condition requires the truncation error 
tends to zero upon (At/ Ax) 2 approach zero as At and Ax approach zero. For this reason, 
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and a much smaller time step than allowed by the above stability limit is implied. This 
concludes that each of the finite difference approximations to the 2-D coupled Burgers' 
equation, the fourth-order explicit Du Fort Frankel scheme and the compact scheme, 
satisfies the consistency condition. Then, the stability of the scheme will be the necessary 
and sufficient condition for convergence, which is true for linear PDE's. But, for the 
present nonlinear PDE's (1), the results of the test cases will verify the convergence, but 
with higher restricted stability limit. 



3 Numerical experiments 

For small value of v, Burgers' equation behaves merely as hyperbolic partial differential 
equation, and the problem becomes very difficult tosolve as steep shock-like wave fronts 
developed, as reported by Kreiss [57]. Therefore, the present higher-order schemes are 
applied to solve problems that are dominated by moderate to severe internal and boundary 
gradients. 



Problem case — 1 



The first test case is the solution of 2-D unsteady coupled Burgers' equations (1) in the 
domain {— 1 < x < 1, < y < ir/6k} with initial and Dirichlet boundary conditions given 
by the exact steady-state solutions, that are set to form moderate to severe internal and 
boundary gradients in the domain [48], as shown in Fig. 2a: 



u s (x,y) 



2 <Pi(x,y) 
Re 4>(x,y) ' 



v s (x,y) 



2 fa(x,y) 
Re 4>(x,y) ' 



where 



4>(x,y) = a + a\x + 
<t>i{x,y) = ai + k 



,fc(x-l) 



+ e 



-k(x-l) 



cos ky, 



e fc(x-l) _ e -k(x-l) 



i(x,y) 



-k 



e fc(x-l) + e -fc(x-l) 



cos ky, 
sin ky. 



(3.1) 



The values of the parameters (ao, a±, k, Re) determine the type of the gradient in the 
computed solutions. We consider three cases in this problem: 

case la. moderate internal gradient with ao = a\ = 110.13, k = 5 and Re = 10; 
case lb. severe internal gradient with ao = a\ = 1.2962 x 10 13 , k = 25 and Re = 50; 
case lc. severe boundary gradient with ao = ci\ = 0.011013, k = 5 and Re = 10. 

The numerical steady-state solutions of the equations (1) have been obtained at time 
= 0.1 for the above three cases, using the present schemes for different grid sizes. Fig. 2b 
shows the computed values of u- velocity component for two different grids; (10 x 5), 
(40 x 20), using the compact ADI scheme. The fourth-order accurate compact ADI scheme 
is capable of producing convergent and stable steady-state solutions with severe gradient 
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even on relatively coarse grid size and large time step size (At = 0.01), in comparison with 
the fourth-order Du Fort Frankel scheme that required finer grid (80 x 40), and smaller 
time step size (At = 0.0002), as shown in Fig. 3. Moreover, the the fourth-order Du 
Fort Frankel scheme exhibits overshoots at the steep gradients especially with a coarse 
grid. This indicates that the fourth-order Du Fort Frankel scheme is unstable, and a much 
smaller time step and grid step sizes than allowed by the linear stability condition are 
required. To test the effect of the initial conditions on the performance of the present 
schemes, initial conditions different from the ones given by equation (3.1), Uij = 1 &; 
Vij = u/um, are used. The results obtained by the compact ADI scheme are stable and 
convergent, while the fourth-order Du Fort Frankel scheme suffers from overshoots at the 
gradients, as shown in Fig. 3b. The computational efficiency of the two schemes has been 
tested by measuring the execution times, using PC-80486 DX2/66, necessary to obtain 
steady-state solutions at time = 0.1, and by computing the two-dimensional error norms 
(E U ,E V ) defined by: 

j N M N M 

1=1 j=l 1=1 j=l 

where (uij,v,ij) represent the numerical computed solutions and (u s ,v s ) the exact solutions. 
Table 1 shows a comparison of the execution times and the error norms of the present 
schemes; the compact ADI scheme and the fourth-order Du Fort Frankel scheme. The 
compact ADI scheme computations with Ai max = 0.01 required about 12.5 seconds, in 
comparison with 26 seconds needed by the fourth-order Du Fort Frankel scheme with 
At max = 0.0002, for the same grid (40 x 20). Moreover, the compact ADI scheme produces 
more accurate and hence more efficient solution, even on a courase grid, as shown in 
Table 1 . This concludes that the fourth-order accurate compact ADI scheme is twice more 
economical and more accurate than the fourth-order accurate Du Fort Frankel scheme. 

Problem case— 2 

In this test case, we consider the solution of 2-D unsteady coupled Burgers' equations (1), 



that is dominated by internal gradients, in the domain {0 < x < 1, < y < 1, t > 0}, 
with the following initial and Dirichlet boundary conditions [50]: 

u(x, y, 0) = sin(7ra;) sin(7ry), (3-3) 

v(x, y, 0) = {sin(7rx) + sin(27rx)}{sin(7ry) + sin(27ry)}, (3.4) 

u(0, y, t) = u(l, y, t) = u(x, 0, t) = u(x, 1, t) = 0, (3-5) 

v(0, y, t) = v(l, y, t) = v(x, 0, t) = v(x, 1, t) = 0. (3.6) 
The computed values of the velocity components u and v, for the case of Re = 1 and 



at different times (t = 0, 0.01, 0.05), on two different grids (40 x 40) , (10 x 10) , using 
the compact ADI scheme, are shown in Fig. 4. Again, the fourth-order compact ADI 
scheme is capable of producing stable and accurate solution with internal gradient on a 
coarse grid (10 x 10), and with comparable accuracy to the solution on a refined grid 
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(40 x 40) . Tables 2-3 show a comparison of the computed solutions using the present two 
schemes, with the previous numerical results of Arminjon and Beauchamp [50], using the 
finite element and the method of lines. The comparison shows that the solutions obtained 
by the compact ADI scheme, and by both the finite element and the method of lines 
coincide to three significant digits in most cases. The behaviour of the solution for u and 
v components at larger time is also considered. The computed solution decreases very 
rapidly to zero at time = 0.5, which is the same result obtained previously by Arminjon 
and Beauchamp [50]. Moreover, the present compact ADI scheme has the advantage over 
the other schemes that it can use a coarse grid and time step size at least twice as larger as 
for the other schemes to get convergent and accurate solution. Fig. 5 shows a comparison 
of the computed solutions for v-velocity at time = 0.01, using the compact ADI scheme, 
and using the fourth-order Du Fort Frankel scheme that required finer grid (40 x 40) and 
smaller time step size (At = 10~ 6 ) to produce stable solutions. Also, the execution times, 
using PC-80486 DX2/66, necessary to obtain stable solutions at time = 0.1 and Re = 1 for 
the present two schemes, are listed in Table 4. The compact ADI scheme computations, 
with Ai max = 10 -3 , required about 85 seconds, while the corresponding computations 
with the fourth-order Du Fort Frankel scheme, with Ai max = 10 -6 , required 780 seconds, 
for the same grid (40 x 40). This indicates that the present fourth-order accurate compact 
ADI scheme has higher computational efficiency for solving the unsteady coupled Burgers' 
equations with internal gradients. 

4 Conclusion 

In conclusion, the fourth-order accurate two-point compact ADI scheme and the fourth- 
order accurate Du Fort Frankel scheme are used to solve the two-dimensional unsteady 
coupled Burgers' equations, for problems that are dominated by moderate to severe in- 
ternal and boundary gradients. The accuracy and the computational efficiency of the 
present schemes are tested. The compact ADI scheme is found to be stable, efficient, and 
with better resolution of steep gradients in comparison with the other scheme, and with 
the previous numerical results of Arminjon and Beauchamp, using the finite element and 
the method of lines [50] . 
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Table 1. Comparison of execution times and error norms for computed solutions of 2-D 
coupled Burgers' equations at time = 0.1 & Re = 50, for problem case-1. 





4-th compact scheme 


4-th Du Fort Frankel scheme 


grind 
points 


execution time 
(seconds) 


-logE u 


-log 


execution time 
(seconds) 


- log E u 


- log E v 


10 x 5 


At = 0.01 
time = 1 


12.0 


11.4 


At = 0.0002 
time = 2.5 


7.4 


5.4 


20 x 10 


At = 0.01 
time = 3 


13.6 


13.3 


At = 0.0002 
time = 7 


5.4 


4.6 


40 x 20 


At = 0.01 
time = 12.5 


13.8 


13.9 


At = 0.0002 
time = 26 


4.8 


4.2 


80 x 40 


At = 0.01 
time = 42 


14.8 


15.4 


At = 0.0002 
time = 100 


4.5 


4.1 



Table 2. Comparison of the values of u and v computed by the Compact ADI scheme 
and Du Fort Frankel scheme at time = 0.01 k, Re = 1, for problem case-2. 



points 


compact ADI scheme 


4-th Du Fort Frankel scheme 




N = 10 
At = 1/800 


N = 20 
1/800 


N = 40 
1/1000 


N = 10 
At = 10" 6 


N = 20 
10~ 6 


N = 40 
10~ 6 


(0.1,0.1) 
(0.2,0.8) 
(0.4,0.4) 
(0.7,0.1) 
(0.9,0.9) 


the velocity u 

0.07320 0.07275 0.07273 
0.27800 0.27803 0.27800 
0.72292 0.72290 0.72285 
0.20542 0.20506 0.20497 
0.07968 0.07955 0.07953 


the velocity u 

0.09549 0.08446 0.07729 
0.30231 0.28314 0.27929 
0.72942 0.72215 0.72183 
0.25000 0.21957 0.20978 
0.09549 0.08785 0.08276 


(0.1,0.1) 
(0.2,0.8) 
(0.4,0.4) 
(0.7,0.1) 
(0.9,0.9) 


the velocity v 

0.43599 0.43662 0.43448 
-0.13444 -0.13131 -0.13148 
1.65503 1.65869 1.65917 
0.06486 0.06337 0.06417 
0.01427 0.01512 0.01476 


the velocity v 

0.80425 0.63005 0.51128 
-0.29180 -0.16784 -0.13866 

1.66391 1.65563 1.65177 
-0.12738 0.00655 0.04585 

0.07771 0.04794 0.02743 
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Table 3. Comparison of the values of u and v computed by the method of lines and the 
finite element method [50] at time = 0.01 k, Re = 1, for problem case-2. 



points 


method of lines 


finite element method 




TV = 10 
At = 1/800 


N = 20 
1/2000 


N = 40 
1/4000 


N = W 
At = 1/2000 


N = 20 
1/1500 


N = 40 
1/2000 


(0.1,0.1) 
(0.2,0.8) 
(0.4,0.4) 
(0.7,0.1) 
(0.9,0.9) 


0.07277 0.07257 0.07253 
0.28887 0.28846 0.28836 
0.72315 0.72205 0.72178 
0.20139 0.20112 0.20106 
0.07956 0.07948 0.07947 


f hp \r£±\r\r t it'\r if 

0.07279 0.07257 0.07252 
0.28867 0.28842 0.28835 
0.72370 0.72210 0.72179 
0.20157 0.20117 0.20107 
0.07951 0.07947 0.07946 


(0.1,0.1) 
(0.2,0.8) 
(0.4,0.4) 
(0.7,0.1) 
(0.9,0.9) 


the velocity v 

0.43857 0.43302 0.43173 
-0.13200 -0.12387 -0.12184 
1.66509 1.65571 1.65335 
0.06137 0.06571 0.06679 
0.01459 0.01372 0.01349 


the velocity v 

0.44130 0.443357 0.43178 
-0.13172 -0.12366 -0.12180 
1.66212 1.65499 1.65316 
0.06306 0.06621 0.06692 
0.01459 0.01367 0.01349 



Table 4. Comparison of execution times for computed solutions of 2-D coupled Burgers' 
equations at time = 0.01 Sz Re = 1.0 for problem case-2. 



grid points 


4-th compact ADI scheme 


4-th Du Fort Frankel Scheme 


10 x 10 


At = 1.25 x 10~ 3 
time = 4 sec 


At = 10~ 6 
time = 35 sec 


20 x 20 


At = 1.25 x 10" 3 
time = 14 sec 


At = 10~ 6 
time = 240 sec 


40 x 40 


At = 10- 3 
time = 85 sec 


At = 10~ 6 
time = 780 sec 
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Figure Captions 

Figure 1: The couputed amplification factors (2.22) of the numerical solution of coupled 
Burgers' equations (1), using the explicit fourth-order accurate Du Fort Frankel scheme 
for different values of |c| and d. 
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Figure 2a: The excact steady state solutions for u and v velocities of the coupled unsteady 
Burgers' equations for problem case-1. 
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Figure 2b: The computed steady state solutions for w-velocity of the coupled Burgers' 
equations for problem case-1, for two different grid sizes; (40 x 20), (10 x 5), using the 
present compact ADI scheme. 
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Figure 3: The computed steady state solutions for u and v velocities of the coupled 
unsteady Burgers' equations for problem case-lb, with internal severe gradient, using 
compact ADI scheme, and Du Fort Frankel scheme. 
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Figure 4: The computed solutions for u and v velocities of the coupled unsteady Burgers' 
equations for problem case-2 at t = 0, 0.01, 0.05, & Re = 1, for two different grid sizes 
(40 x 40), (10 x 10), using the present compact ADI scheme. 
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Figure 5: Comparison of the computed solutions of f-velocity component of the coupled 
unsteady Burgers' equations for problem case-2, at time = 0.01 k, Re = 1, using the 
present compact ADI scheme, and Du Fort Frankel scheme. 
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